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Abstract. The heliosphere is predominantly a partially ionized plasma that con- 
sists of electrons, ions and significant neutral atoms. Nonlinear interactions amongst 
these species take place through direct collision or charge exchange processes. These 
interactions modify linear and non linear properties of the plasma. In this work, we 
develop a one dimensional linear theory to investigate linear instabilities in such 
system. In our model, the electrons and ions are described by a single fluid com- 
pressible magnetohydrodynamic (MHD) model and are coupled self-consistently to 
the neutral fluid via compressible hydrodynamic equations. The coupling is me- 
diated by the charge exchange process. Based on our self-consistent analysis, we 
find that the charge exchange coupling is more effective at larger length scales, and 
Alfven waves are not affected by the charge exchange coupling. By contrast, the 
fast and slow waves are driven linearly unstable. 



1. Introduction 

Magnetohydrodynamic (MHD) plasma admits fundamental modes in the form of 
Alfven, fast and slow waves. These modes are ubiquitously observed in many lab- 
oratory [5], space [6, 14] and astrophysical plasmas [1]. The Alfven, fast and slow 
waves are essentially electromagnetic oscillations in a plasma that propagate along 
an ambient or guide magnetic field. Their propagation characteristic depends on 
the orientation of their wave vector relative to the velocity and magnetic field fluc- 
tuations. For instance, wave vector associated with Alfven waves is orthogonal to 
the velocity field fluctuations whereas the fast and slow waves have the parallel 
wave vector relative to the velocity field fluctuations [14, 15]. 

The propagation properties of the Alfven, fast and slow waves have long been 
known and well studied in the MHD plasma (in isolation). But their evolution in the 
presence of complex interactions or other species, such as neutral atoms/particles, 
is far less explored because of untractable nature of analytic as well as numerics. 
In the context of partially ionized (consisting of ions, electrons and significant neu- 
trals) space and astrophysical plasmas, these waves interact with the neutral gas 
and govern numerous properties. Kulsrud & Pearce (1969) in the context of cosmic 
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ray propagation noted that the interaction of a neutral gas and plasma can damp 
Alfven waves. Neutrals interacting with plasma via a relative drag process results 
in ambipolar diffusion [Oishi & Mac Low 2006] . Ambipolar diffusion plays a crucial 
role in the dynamical evolution of the near solar atmosphere, interstellar medium, 
and molecular clouds and star formation. Oishi & Mac Low (2006) investigated 
that the ambipolar diffusion sets a characteristic mass scale in molecular clouds. 
They found presence of structures below the ambipolar diffusion scale because of 
the propagation of compressive slow mode MHD waves at smaller scales. Leake et 
al (2005) showed that the lower chromosphere contains neutral atoms, the existence 
of which greatly increases the efficiency of wave damping due to coUisional friction 
momentum transfer. They noted that Alfven waves with frequencies above 0.6Hz 
are completely damped and frequencies below 0.01 Hz are unaffected. They under- 
took a quantitative comparative study of the efficiency of the role of (ion-neutral) 
coUisional friction, viscous and thermal conductivity mechanisms in damping MHD 
waves in different parts of the solar atmosphere. It was pointed out by the authors 
that a correct description of MHD wave damping requires the consideration of all 
energy dissipation mechanisms through the inclusion of the appropriate terms in 
the generalized Ohms law, the momentum, energy and induction equations. Padoan 
et al (2000) calculated frictional heating by ion-neutral (or ambipolar) drift in tur- 
bulent magnetized molecular clouds and showed that the ambipolar heating rate 
per unit volume depends on field strength for constant rms Mach number of the 
flow, and on the Alfvenic Mach number. 

In the context of molecular clouds, the role of ion-neutral collision has been inves- 
tigated by Balsara (1996). It was found that momentum of the plasma is governed 
predominantly by the slow comoving neutrals that tend to dissipate Cloud's mag- 
netic field. Balsara (1996) described a linear analysis by considering neutral momen- 
tum in the presence of gravitational field of the Clouds that dominates the plasma 
momenta because of their higher masses. Hence they ignored the inertial terms in 
plasma momentum equation. Balsara (1996) reported that in super Alfvenic regime, 
the slow waves propagate without significant damping near short wavelengths, while 
the fast and Alfven waves undergo rapid damping. 

In addition to the neutral gas interacting collisionally with the plasma, it under- 
goes charge exchange process in the heliosphere and in the local interstellar medium 
(ISM) [17, 12]. Various regions of interactions in the heliosphere are depicted in Fig 
(1). The charge exchange process knocks off electron from a neutral and makes it 
a ion. This electron is captured by the plasma ion which then turns into a neutral 
atom. This process conserves the density and but not momentum and energy of 
the plasma protons and H neutral gas. In the local interstellar medium (LISM), 
just beyond the heliosphere, the low density plasma and neutral hydrogen (H) gas 
are coupled primarily through the process of charge exchange. Thus on sufficiently 
large temporal and spatial scales, a partially ionized plasma is typically regarded 
as equilibrated. Nonetheless it remains to be seen as to what happens to these 
MHD modes in the presence of charge exchange process, both in the inner and 
outer helioshpere regions. We address this question in this paper based on a linear 
theory that describes the coupled plasma-neutral system mediated by the charge 
exchange process within the realm of heliosphere. While obviously of interest to the 
solar wind-LISM interaction [17, 12], this work will have interesting implications 
for any shock, structures, or turbulent fluctuations embedded in a partially ionized 
medium, ranging from supernova shocks (SNR) to cometary bow shocks. 
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Figure 1. Schematic overview of different regions in the global heliosphere. The solar 
wind emanating from the Sun propagates outward and interacts with partially ionized 
interstellar gas predominantly via charge exchange, and creates pick up ions (PUIs). At 
the termination shock (TS), the supersonic SW decelerated, heated, compressed becoming 
subsonic, in the heliosheath, and again interacts with interstellar neutrals via charge ex- 
change before it reaches heliopause (HP). The subsonic SW flows down into the heliotail. 
During its journey from the Sun to the IfP, the solar wind plasma develops multitude of 
length and time scales that interact with the partially ionized interstellar gas, TS, and 
nonlinear structures develop in a complex manner. 

In section 2, we describe the governing equations of the coupled plasma- neutral 
system in the heliosphere plasma. Importance of the charge exchange is discussed. 
Section 3 deals with the linear theory of MHD modes in the presence of charge 
exchange. Here the plasma and neutral fluids are coupled via charge exchange. 
Because of the analytic complexities, we focus only on one dimensional calculation. 
The solution of linear dispersion relation for the coupled plasma-neutral system is 
described in section 4. In section 5, we summarize our results. 



2. Model equations 

Within the paradigm of the heliosphere plasma, the plasma and neutral fluid remain 
close to thermal equilibirium and behave as Maxwellian fluids. They are coupled 
through charge exchange forces. There may additionally be direct collision between 
the two components. We however ignore it for our current calculation. In general, 
the fluid model describing nonlinear turbulent processes, in the presence of charge 
exchange and collision, can be described in terms of the plasma density (p) , velocity 
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(U), magnetic field (B), pressure (P) components according to the conservative 
form 



dF 



dt 



r + V • = Qp,„, 



(2.1) 
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Note the presence of the source terms Q that couple the plasma sclf-consistcntly to 
the neutral gas. The above set of plasma equations is supplimented by V • B = and 
is coupled self-consistently to the neutral density (A'"), velocity (V) and pressure 
(P„) through a set of hydrodynamic fluid equations, 
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Equations (2.1) to (2.2) form an entirely self-consistent description of the coupled 
plasma-neutral turbulent fluid in a partially ionized medium. 

Several points arc worth noting. The charge-exchange momentum sources in the 
plasma and the neutral fluids, i.e. Eqs. (2.1) & (2.2), arc described respectively by 
terms QM,p,n(U, V, p, N, T,^, Tp) and QM,n,p(V, U, p, N, Tn,Tp). These expressions 
arc described in Pauls ct al (1995) and Shaikh & Zank (2008, 2010). A swapping of 
the plasma and the neutral fluid velocities in this representation corresponds, for 
instance, to momentum changes (i.e. gain or loss) in the plasma fluid as a result 
of charge exchange with the neutral atoms (i.e. QM.p.ri(U, V, p, A'^, T„, Tp) in Eq. 
(2.1)). Similarly, momentum change in the neutral fluid by virtue of charge exchange 
with the plasma ions is described by QM,n,p(V, U, p, N, T„, Tp) in Eq. (2.2). 

The underlying coupled fluid model can be non-dimcnsionalizcd straightforwardly 
using a typical scale-length (£o), density (po) and velocity {vq). The normalized 
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Figure 2. Solution of linear dispersion relation in the smaller k (large scale) regime char- 
acterizes different MHD and hydrodynamic modes in the coupled plasma-neutral system. 
Upper panel describes real part of the frequency where forward/backward fast and slow 
waves are shown. Growth and damping of linear instability are shown in the lower panel, 
where 7 = root corresponds to the sound wave. The latter is stable at lower k's. 



plasma density, velocity, energy and the magnetic field are respectively; p = p/ po, tJ = 
U/i^o, Pp = Pp/poVQ,^ = B/wqVPo- The corresponding neutral fluid quantities are 
N = N/pojV = V/votPu = Pu/poVq- The momentum and the energy charge- 
exchange terms, in the normalized form, are respectively Q™ = Qm^o/po^O'Qe ~ 
Qe^o/ PqVq ■ The non-dimensional temporal and spatial length-scales are t = tvo /£o,^ — 
x/£o- Note that we have removed bars from the set of normalized coupled model 
equations (2.1) & (2.2). 
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3. Linear theory 

To avoid mathematical complications associated with the coupled plasma neutral 
system, we develop a one dimensional linear theory to investigate growth/damping 
of fundamental AlHD modes in the hcliosphcrc plasma. One dimensional theory 
offers not only a better insight into the complex (charge exchange) interaction 
processes, but it also is easier (compared to two dimensional treatment) to deal with 
analytically. Furthermore, it describes dominant features of the characteristic modes 
of the coupled plasma-neutral system. In our linear theory, the charge exchange 
term Qm in the momentum equation is proportional to vpiU — V) [Florinski et 
al 2005] where i/ is the charge exchange frequency. This is a simpler form of the 
charge exchange force. For more complicated structure of the charge exchange force, 
we employ higher order terms for the nonlinear simulations. This is described in 
Shaikh & Zank (2010) where nonlinear stage of mode coupling instability in two 
dimension was investigated. Our results, to be deduced below, help explain some of 
the features of the nonlinear simulations described in Shaikh & Zank (2010). The 
one-dimensional plasma equations in the normalized form are described as; 
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Here u = aU*N is charge exchange frequency that depends on the charge exchange 
cross section (cr) and the relative difference in the speed of neutral and plasma 
components. The plasma thermal speed is Vth = \J kTi/rrii and 7 is the specific heat 
ratio (it is 5/3 in our model). Note that the sign of the momentum charge exchange 
force is negative. This means that charge exchange process tends to decrease plasma 
momentum relative to the neutral fluid. The magnetic field is unaffected by the 
charge exchange process and it continues to be frozen into the plasma component 
of the flow. Similar equations for the neutral fluid density {N)^ velocity {V) and 
pressure (P„) can be described as follows. 
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(3.6) 



(3.7) 



Unlike plasma, the momentum charge exchange force corresponding to the neutral 
fluid is positive. This leads to an enhancement in the neutral momentum relative 
to the plasma fluid. The rate of decay or enhancement of both the plasma and 
neutral momentum also depends on the relative speed between the two species. It 
is evident from Eqs. (3.1) to (3.7) that the plasma and neutral fluids are coupled 
through charge exchange sources. 
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Figure 3. Linear instability in the moderately higher k regime. Backward and forward 
fast/slow modes are shown respectively by the lower and upper curves in the upper panel. 
Interestingly the modes that decay in the small k regime, tend to grow eventually when k 
is order unity. This is shown in the lower panel. 

We next carry out a linear analysis by assuming fixed (or constant) background 
quantities around which the small amplitude perturbed quantities are expanded. 
The amplitude of the linear perturbation is assumed to be smaller compared to the 
background quantity for a valid description of the linear theory. All the perturbed 
variables are expressed in the form of waves with frequency uj and wavenumber k. 
The plasma quantities are described as 

p{x,t)=po + Pk e'"*, 



P{x,t)^Po + Pk e'''- 
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B{x,t)=Bo + Bk e*'^"-^"*, 

where pk{x,t),Uk{x,t)Pi.{x,t) and Bk(x,t) are respectively small amplitude den- 
sity, velocity, pressure and magnetic field fluctuations. Note that the linear fre- 
quency LJ is essentially a complex quantity i.e. a; = cj^ + ^7 where the linear part 
corresponds to the wave. By contrast the imaginary part describes growth (if pos- 
itive) or damping (if negative) associated with the linear perturbations. We are 
interested in determining 7(fc) in our work. 

The corresponding linear expansion for the neutral fluid is 

N{x,t) = NQ + Nk 6*'=^-'"*, 
V{x,t) = Vo + Vk (^^^-'^\ 

where A'^g, Vo,P„„ are respectively the background constant neutral density, veloc- 
ity and pressure, and Nk,Vk,Pnk are the amplitude of corresponding perturbed 
variables. 

On linearizing Eqs. (3.1) to (3.7), we obtain 

-iojpk + ik{poUk + UoPk) = 0, (3.8) 
-iujpoUk + poUoikUk + ikPk = -y{poUk + PkUo - PoVk - PfeVb), (3.9) 

-iioPk + ikUoPk + ikPolUk = upk ^^{Uo-Vo)^ + v^,,-p^ , (3.10) 

L 2 2po_ 

LO = kUo, (3.11) 

-iioNk + ik{NoVk + VoNk) = 0, (3.12) 

-iujNoVk + NoVoikVk + ikPn^ = v{poUk + PkUo - poVk - PkVo), (3.13) 



-iu>Pn„ + ikVoPn^ + ikPnojVk = vpk 



(3.14) 



It is noted that the Alfven wave frequency, oj = kUo, is decoupled trivially from 
the neutral fluid despite charge exchange coupling interactions. This indicates that 
Alfven waves are not affected by the charge exchange interactions. This result can 
be elucidated as follows. Firstly, the charge exchange process tends to swap the 
plasma ions with the neutral particles. During this process, the wave activity is 
not tampered by the charge exchange. Secondly, the magnetic fleld is unperturbed 
by the charge exchange process and it continues to remain frozen in the plasma 
component of the flow. Hence the linear frequency of Alfven waves continues to 
remain unchanged. However, this scenario can change during the nonlinear evolu- 
tion in that Alfven waves are shown to exhibit modulation [Shaikh & Zank 2010] 
mitigated by rapid nonlinear growth and damping that occur alternatively in the 
coupled plasma-neutral system. 

After eliminating all the perturbed variables, we get the fifth order linear disper- 
sion relation with complex coefiicients, 



A5 J' -I- ^4 -h A3 -h A2 -I- ^1 w Ao = 0. 



(3.15) 



Note further that the sixth order term in the above polynomial exists separately 
in the form of the Alfvenic branch (i.e. w = kUo) which is shown to decouple 
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because the magnetic field, critical in exciting Alfven waves, is unaffected largely 
by the charge exchange process. The magnetic field fluctuations are embeded into 
the plasma component and are affected only indirectly by the plasma momentum 
flows. Nonetheless, the complex coefficients in Eq. (3.15) are given below; 

A5 = No, 

A4 = ipoiNou + iy)-2 VoNok - poUoNok - 2 UoNok, 

As = ipoUo{-Noku - 2ku) + ipQVQ{-^Nokv - ku) 

+iplUokiy + ji-NoPok'' - P„ J + Po(2 UoVoNok 
+2 U^Nok"") + V^Nok^ + 4 UoVoNok'' + U^VoNQk\ 



A2 = PoNok^i^6p - pUoV„^Nok^ - mVoNok^ - U^Nnk"" 
+ipoUoVoiSNok^iy + 2k^v) + ipoW^Nok^u + 
ipliUoVok^v - 2Uieu) + ip^oU^e) 
+^i{-poPok'^v - PnoPo) + 2-fVoNoPok^ 
+'yUoNoPok^ + PnoipoUok + 2Uok)'y - 2UoV^Nok^ - 2USVoNok\ 



Ai = ipo{-2VoNok^u5p - V^Nok^v - WoV^Nok^v - 
U^Vok^v) + ipliU^k'^v - 2U^Vok^iy) + pUv'^S^ 
-k^v'^Sp) + 7ipo(VbPofc^J^ + UoP^k^v) 
+iP„,Poi2Uok - lyUok + Vok) - iVi^N^P^k^ - 2iUoVoNoPok^ 
+iPno{-2poUik' - U^k") + po{2U^VoNok^ + 2UivSNQ) 

+U^ViNok^ + 'yPnoPo, 



Ao = -poU^VoNok^ - ^^PnoUoPok 
+'yiPnoPo[-k''iy5p + U^k^iy - 1) - UoVok^] 
-ipoUoVoPoiyk* + jUoVo^NoPok" + PnoPoU^k"" 
+ipoV^k^i^6p + ipoUoViNak'^v + voC/q KtfcV 
+pl{Vok^v^5p - Uok6„iy^), 



6p 



, ■ -,r \2 , 2 -^0 

— (Uc-V„) +«.»-2^ 



4. Solution of the dispersion relation 

We next describe the solution of linear dispersion relation, Eq. (3.15), for the cou- 
pled plasma-neutral system in the heliosphere. To access a broader spectrum, we 
solve Eq. (3.15) in three distinct regime, namely the short length scale {k/kc > 1), 
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Figure 4. Larger modes k/kc > 1 (smaller scales) show predominant propagation of 
forward fast /slow waves as shown in the upper panel. Smaller scales show a decoupling 
with the sound waves. Lower panel shows unstable modes in the coupled plasma-neutral 
system for larger modes. Sound waves are found to be stable for the smaller scales. 

longer length scale {k/kc < 1) and the scales comparable to the charge exchange 
length (k/kc ~ 1), where kc = y/cmo = 2t:/Ic is the charge exchange mode (Ic charge 
scale length) . Charge exchange interactions not only modify the propagating char- 
acteristic of the coupled plasma-neutral modes, but they also alter the growth and 
damping rates. 

We first examine the affect of charge exchange coupling on the longer (than charge 
exchange length Ic) length scale part of the spectrum. The real part of the complex 
frequency leads to the sound wave, forward and backward propagating fast/slow 
magnetosonic waves. This is shown in Fig (2). Alfvenic branch is not shown in 
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this figure. Note that the sound waves in the neutral fluid coexist with the plasma 
modes due to the charge exchange coupling. Furthermore, the sound waves are not 
destablized by the charge exchange at lower fc's. Interestingly, the charge exchange 
interactions drive large scales unstable. It is evident from Fig (2), lower panel, that 
the linearly unstable damping modes coexist with the growing (or unstable) modes. 
By virtue of this, the mode coupling interactions turn out to be far more complex. 
An implication of this instability is to lead to an alternate growth and damping of 
the fast/slow waves. This instability was thought to govern nonlinear modulation 
of the fast/slow wave in the coupled plasma-neutral system reported in Rcf [15] by 
us. It is further noted that the sound waves begin to decouple from the MHD modes 
for higher fc > 1 modes. This is an indication that the small scales are not coupled 
efficiently by the; charge exchange interactions. This is further consistent with our 
nonlinear two dimensional simulations [14] that show a weak coupling between the 
small scale turbulent fluctuations and charge exchange scales. 

We next access a moderately higher wavenumber regime in the spectrum. This 
regime, depicted in Fig (3), exhibits an entirely different characteristic compared 
to the one described in Fig (2). Clearly, the forward propagating fast and slow 
magnetosonic modes dominate over the slow (below zero roots) ones in the higher 
k spectrum where k > kc (but not k kc). By contrast, the unstable modes 
(lower panel in Fig 3) exhibits interesting feature. In that, the linearly unstable 
fast/slow modes (i.e the second curve from the top) shows initial growth (i.c 7 > 0) 
for k < 1.25. This characteristic behavior reverses for k > 1.25. Converse is true 
for the third (from the top) curve. On the other hand, the top and lowest curves 
corresponding to the forward and backward propagating fast/slow modes show 
respectively the growth and damping of the instability. 

Finally, to investigate the asymptotic behavior of these linearly stable/unstable 
modes, we extend our analysis to span the higher k modes in the spectrum. This 
is shown in Fig (4). It becomes clear from this figure that the linear wave propaga- 
tion (upper panel) show a considerable shift in the frequency towards the forward 
propagating branch {wr > 0), whereas the frequency associated with the back- 
ward propagation is found to be reduced considerably. Interestingly, the fast/slow 
branches associated with the magnetosonic waves appear to have merged for the 
smaller scales in the fc-spectrum. The unstable modes, (Fig 4, lower panel), on 
the other hand exhibit a predictable behavior. It is found that the forward prop- 
agating fast/slow modes are driven unstable by the charge exchange coupling in 
the coupled plasma-neutral system. On the other hand, the backward propagating 
fast /slow modes are damped by the charge exchange interactions. This is described 
by the lower panel in Fig (4). It is further observed that the sound waves, in the 
higher k regime, are decoupled from the predominant MHD modes. This means 
that charge exchange do not mix up these modes at relatively small characteristic 
length scales. 

5. Summary and conclusion 

A major outcome from our analysis is that charge exchange interactions in the par- 
tially ionized gases, dominated by the coupled plasma-neutral system, in the he- 
liosphere do not directly modify the propagation characteristic of the linear Alfven 
waves. This does not necessarily mean that their nonlinear evolution is also unaf- 
fected. As a matter of fact, the nonlinear evolution of Alfven waves is influenced 
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by their mode coupling interaction with the fast/slow magnetosonic waves. This 
is discussed by Shaikh & Zank (2010). We however find that the charge exchange 
interaction influence fast/slow mode by hnearly driving them unstable. 

The linear analysis described in the paper is extremely useful in explaining some 
of the features of our nonlinear simulations [14, 15]. The latter show that charge 
exchange modes modify the helioshperic turbulence cascades dramatically by en- 
hancing nonlinear interaction time-scales on large scales. By contrast, small scale 
turbulent fluctuations are unaff'ected by charge exchange modes which evolve like 
the uncoupled system as the latter becomes less important near the larger k part 
of the helioshperic turbulent spectrum. This tends to modify the characteristics of 
helioshperic plasma turbulence which can be significantly different from the Kol- 
mogorov phenomenology of fully developed turbulence. 

The onset of modulational instability in Ref [15] is triggered essentially by lin- 
ear instability process. Our linear stability analysis indicates that the underlying 
coupled plasmaneutral system possesses fast/slow unstable modes. These modes 
account for the growing as well as damping of Alfvnic and fast/slow compressive 
waves. The latter leads to a concurrent damping of the fast/slow compressive mode 
preceding the growth of Alfvenic mode occurs. When the two modes reach their 
extremal rise or fall, they reverse their behavior. The fast/slow compressive mode 
begins to rise at the expense of damping of Alfvenic mode. This process continues 
nonperiodically and repeats itself in time thereby exhibiting a nearly oscillatory 
behavior. It then follows that the two modes, i.e. Alfvenic and fast/slow compres- 
sive, regulate each other in a predatorprey manner. Dynamically, the Alfvenic mode 
grows at the expense of fast/slow mode. When Alfvenic mode reaches its maximum 
amplitude, where the fast/slow mode remains at its lowest magnitude, the latter 
eats up the Alfvenic mode and vice versa. This process is repetitive. 

Our results should find application to a variety of astrophysical environments in 
which a partially ionized plasma is typical. Examples include the outer heliosheath 
formed by the interactions of the solar wind with the local interstellar medium, the 
magnetic collapse of molecular clouds and star formation and the general transfer 
of energy in partially ionized plasmas surrounding other stellar systems [16]. 
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